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Stabilization of a linear Korteweg-de Vries equation 
with a saturated internal control 


Swann Marx^, Eduardo Cerpa^, Christophe Prieur^ and Vincent Andrieu^ 


Abstract —This article deals with the design of saturated 
controls in the context of partial differential equations. It is 
focused on a linear Korteweg-de Vries equation, which is a 
mathematical model of waves on shallow water surfaces. In this 
article, we close the loop with a saturating input that renders 
the equation nonlinear. The well-posedness is proven thanks to 
the nonlinear semigroup theory. The proof of the asymptotic 
stability of the closed-loop system uses a Lyapunov function. 

I. Introduction 

In recent decades, a great effort has been made to take 
into account input saturations in control designs (see e.g [24] 
or [10]). Indeed, in most of systems, actuators are limited 
due to some physical constraints and the control input has 
to be bounded. Neglecting the amplitude actuator limitation 
can be source of undesirable and catastrophic behaviors for 
the closed-loop system. The standard method follows a two 
steps design. First the design is carried out without taking 
into account the saturation. In a second step, a nonlinear 
analysis of the closed loop system is made when adding 
the saturation. In this way, we often get local stabilization 
results. Tackling this particular nonlinearity in the case of 
finite dimensional systems is already a difficult problem. 
However, nowadays, numerous techniques are now available 
(see e.g. [24], [25], [23]) and such systems can be analyzed 
with an appropriate Lyapunov function and a sector condition 
of the saturation map, as introduced in [24]. 

To the best of our knowledge, there are few papers 
studying this topic in the infinite dimensional case. Among 
them, we find [12] and more recently [18], where a wave 
equation equipped with a saturated distributed actuator is 
considered. Note that saturation function can be defined 
with a sign function, which is also used in sliding mode 
control design theory. The interest reader can refer to [17], 
[8], where a wave and a reaction-diffusion equations are 
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stabilized with a sliding mode controller. The present paper 
aims at contributing to the study of the saturated input case 
in the framework of partial differential equations. 

Let us note that in [22] the case of a priori bounded 
feedback is studied for abstract linear systems. To be more 
specific, for compact control operators, some conditions 
are derived to deduce, from the asymptotic stability of 
an infinite-dimensional linear system in abstract form, the 
asymptotic stability when closing the loop with saturating 
controller (see [22, Theorem 5.1] for a precise statement of 
this result). The aim of our article is to study a particular 
partial differential equation without seeing it as an abstract 
control system and without checking the very specific as¬ 
sumptions of [22]. 

The Korteweg-de Vries equation (KdV for short) 

yt + Vx + yxxx + yyx = ^, (l) 

is a mathematical model of waves on shallow water surfaces. 
Its controllability and stabilizability properties have been 
deeply studied in the case with no constraints on the control, 
as explained in [3], [6], [20]. In this article, we focus on the 
following controlled linear KdV equation 

yt + yx + yxxx + / = 0, (a:, t) e [O, L] x [O, -boo), 
y{i^ 0) = y{t, L) = y^(t, L) = 0, tG [0, -bcxD), (2) 
2/(0,a;) = yo{x), x G [0, L], 

where y stands for the state and / for the control. As studied 
in [19], if / = 0 and 

then, there exist solutions of (|^ for which the energy does 
not decay to zero. For instance, if L = 27r and yo = 1 — 
cos{x) for all x G [0,L], then y{t,x) = 1 — cos(a;) is a 
stationary solution of conserving the energy for any time 
t. In the literature there are some methods stabilizing the 
KdV equation (|^ with boundary [5], [4], [13] or internal 
controls [16], [15]. Here we focus on the internal control 
case. In fact, as proven in [16], [15], the feedback control 
f{t,x) = a{x)y{t,x), where a = a{x) is a positive function 
whose support is a nonempty open subset of (0,L), makes 
the origin an exponentially stable state. 

The question we want to address is the following. Given a 
feedback control / stabilizing the equation, what do we get 
if we saturate it? Is the equation still stable? We deal with 
the case in which f{t, x) = ay{t, x) with a positive constant 




a and we show that the origin is asymptotically stable for 
the closed-loop system with a saturated input. 

This article is organized as follows. In Section |I^ we 
present our main results about the well posedness and the 
stability of this equation in presence of saturation. Section 
111 is devoted to prove these results by using the nonlinear 


semigroup theory and Lyapunov techniques. In Section IV 


we give some simulations of the equation looped by a satu¬ 
rated feedback. Section |V] collects some concluding remarks 
and possible further research lines. 

Notation: yt (resp. y^;) stands for the partial derivative of 
the function y with respect to t (resp. x) (this is a shortcut for 
resp. ^). (resp. 3) denotes the real (resp. imaginary) 
part of a complex number. Given L > 0, \\ ■ ||l2(o_l) 
denotes the norm in L^{0,L) and H^{0,L) is the set of all 
functions u G L^(0, L) such that Ux,u. 


XX •) ^XXX 


G L2(0,L). 

Finally Hq{0,L) is the closure in L^(0, L) of the set of 
smooth functions that are vanishing at cc = 0 and x = L. 
It is equipped with the norm \MHi( 0 ,L) ■= fo \Ux\‘^dx. 
The associate inner products are denoted {■,-)l^(o,l) 
(•r)ffi(o.i)' denotes the set iT£(0,L) := {w S 

r(;(0) = w{L) = w'{L) = 0}. 


II. Main results 

For any a > 0, if we take f{t,x) := ay{t,x) in (|^, then 
we get that the equation is stabilized. Indeed, any solution 
of 

yt + Vxxx +yx+ay = Q, 

yit,o) = y{t,L) = 0, (3) 

yx (f 1 -^) d; 


satisfies 

1 d 

2 dt jQ 


\y(t,x)\'^dx = -^\y^(t,0)\'^-a \y{t,x)\'^dx 

<-af \y{t,x)\'^dx, 

Jo 


(4) 


We can write the KdV equation controlled by a saturated 
control as follows 

{ yt + yxxx +yx+ asat(y) = 0, 
y{t,0) =y{t,L) =0, 
yx{t,L) = 0, 

2/(0, a;) = yoix). 

Let us state the main results of this paper. 

Theorem 1 (Well-posedness). For any initial condition yo G 
iL£(0,L), there exists a unique strong continuous solution 
y : [0,(30) ^iTl(0,L) to § that is continuous from [0, oo) 
to and continuously differentiable from [0,(X)) to 

L2(0,L). 

Moreover, for any initial condition yo in L^(0,L), there 
exists a unique weak solution y : [0,(X)) -G L^{0,L) to Q 
that is continuous from [0,(X)) to L^(0, L). 

Theorem 2 (Asymptotic stability). For any constant a > 0, 
the equation 0 is globally asymptotically stable. More pre¬ 
cisely, the following property holds. For any initial condition 
yo in L^{0,L), the weak solution to 0 satisfies 

\\y{t, ■)\\l^{o,l) < ||2 /o||l 2 (o,l), Vf > 0, (8) 

together with the attractivity property 

\\y{t,-)\\L^{o,L)-t 0, ast-i’OO. (9) 

Remark 1. The exponential stability of the closed-loop 
system with a saturating control is an open problem for the 
KdV equation. 

III. Proof OF Theorems[I]and[2 
A. Well-posedness (Theorem 0 
Let A denote the operator 

Aw = {—w' — w"' — asat(w)) 

on the domain D{A) C L^(0,L) dehned such that DI^A) := 
Hli0,L). 

Lemma 1. Operator A is closed. 

Proof. Let {u„}„gN be a sequence in D{A) such that 


which ensures an exponential stability with respect to the 
L^(0, L)-norm. Note that the decay rate can be selected as 
large as we want by tuning the parameter a. Such a result is 
refered to as a rapid stabilization result. 

Let us assume now that the control is constrained and that 
we have to consider the following feedback law 


lim Un = u (10) 

n—f-\-oo 

and 

lim Aun = V (11) 


for some u,v G L^(0, L). To prove that A is closed, we have 
to prove that u G D(A) and that Au = v. Let us note that 


f{t,x) = a-sat{y{t,x)) (5) 

where the function sat is dehned by 

ttniin if S ^ rijnin, 

5 if t^min ^ ^ ^ Umax, (b) 
Umax if S > rtniax- 

To ease the lecture, we assume same levels of saturation, 
which means that Umax = Wmin = Mo- 


A:wG D{A) C L2(0, L) ^ {-w' - w'") G L) 

( 12 ) 

is already closed in D{A). Moreover, we know that the 
function sat is globally Lipschit 2 j^ Thu^A is closed. □ 

’indeed, we know from [11, Page 73] that for all {y,y) S L'^{0,Lfi 
and for all x G [0, L], |sat(r/(a;)) — sat(^(fc))| < \y(x) — y(x)\. Thus 
we get ||sat(t/) - sat(y)||^2(o,L) < Wv “ vWL^io.L)- 

^Given Ti closed and T 2 globally Lipschitz, and u„ u and (Ti -t 
T 2 )u„ -a w, we have \TiUn -t T 2 U — ml < |TiUn -t T 2 Un — ui| -t 
\T 2 Un — T 2 u\. Thus, the left member of the inequality is bounded by a 
term which converges to 0. 




Lemma 2. Operator A is dissipative. 

Proof. Let us consider 

satc(s) := sat(iH(s)) +jsat(rj(s)) (13) 

which we will denote by sat(s) to ease the notation. 
Given u,u G D{A), we have that 


ip{u, u) := {Au - Au,u- u)l^o,l) 


is equal to 


ip{u,u)=— / {u"\x) + u'(x) + asat{u{x)) 

Jo 


— {u"'{x) + u'{x) + asat(M(a:))) .{{u — u){x))dx 

_ 

= — a {s at (u) — sat {u)){u — u)dx 

Jo 

— j {u'" — u"'){u — u)dx 

Jo 


(14) 


u'(0)|2 < 0. (15) 


(16) 


Integrating by parts (u'" — u"')(u — u), we get 

=- 

Then we have 

IH {{Au - Au,u- u)l^o,l) } < 

— aDl I J {sat{u) — sat{u)){u — . 

By definition of the saturation function, we get that for all 

(s, s) G 

93 {(sat(s) — sat(s))(s — s)} > 0. (17) 

Thus, thanks to the positivity of a, we get that 

‘^{{Au-Au,u- (18) 

which means that the operator A is dissipative. It concludes 
the proof of Lemma □ 

In order to conclude the proof of the well-posedness, we 
have to verify whether the operator A generates a semigroup 
of contractions which will be denoted in the following by 
S{t). Following [14], we see that it is enough to prove that 
for all A > 0 sufficiently small 


D{A) C Ran(/ - \A), 


(19) 


where Ran stands for the range and I for the identity 
operator. In other words, for each u G D{A), there exists 
u G D{A) such that 


(/ — XA)u = u, 


( 20 ) 


which is equivalent to prove the existence of a solution of a 
nonhomogeneous nonlinear equation in the u-variable with 
boundary conditions as considered in the following lemma. 


Lemma 3. Let us introduce X := j. If a is strictly positive 
and A > 0, then there exists u solution of 


Xu + u" + u + asat{u) = Xu, 
u(0) = u{L) = u'{L) = 0. 


( 21 ) 


Proof. The proof of this lemma follows from classical tech¬ 
nics (see e.g. [14, Page 179]) and uses the Schauder fixed- 
point theorem (see e.g. [6, Theorem B.19,]). 

First, following [7], let us focus on the spectrum of A 
which is defined by (12i. Since the operator A has a compact 
resolvent, its spectrum denoted by a {A) consists only of 
eigenvalues. Futhermore, the spectrum is a discrete subset 
of iM.. 

Since A belongs to IR+, then A ^ o'{A). Hence (^ — /A) 
is invertible and there exists a unique function z = z{x) 
solution of 


\ Xz + z'" + z' = g, 

[ z{0) = z{L) = z'{L) = 0, 

where g{y) := —asat(y) -I- Xu. 

Then we can focus on the map 

T:L^{0,L) -G- L^{Q,L) 
V '—^ = T{y) 


( 22 ) 


(23) 


where z = T{y) is the unique solution to (22i. We define 

C = {uG Hl{0, L)/ ||w|l^i(o,L) < M} (24) 

where M > 0. From the theorem of Rellich (see [2, Theorem 
9.16, p. 285]), the injection of Hq{0,L) in L^{0,L) is 
compact, then C is bounded in Hq{0,L) and is relatively 
compact in Lf{0,L). Moreover, it is a closed subset of 
L^(0, L). Thus 67 is a compact subset of L^(0, L). In 
order to apply the Schauder theorem, we have to prove that 
T{L^{0, L)) C 67 for a suitable choice of M >0. We 


multiply the first line of ( 221 by z and then integrate between 


0 and L. After some integrations by parts, we get 

pL pL 

M\A\l'^(o,l) = - / zz"'dx- / zz'dx 

JO JO 

pL pL 

— a s at (y) zdx + X / uzdx 


nL pL 

sat {y)zdx + X / uzdx 


= - z'{0f - a 

Jo 

pL pL 

< — a s at (y) zdx + X / uzd 
Jo Jo 


(25) 


The Young inequality leads us to the following inequality 


2^11^1112(0,L) <aei||sat(2/)||i2(0,L) + —\\z\\l^o,l) 
+ ^II^IIl2(0,L) + ^£2||M|li2(0,L) 

where ei,e 2 >0 are to be chosen later. 


( 26 ) 












The function sat(-) being bounded, we get 


Since sat is an odd function, thus, for all (f,x) S K+ x 
[0,L] 


X 


, ^ £1 £2 , 


I^IIl2(o,l) < aeiLuQ + A£2||w||l2(o,l) 

We choose ei and £2 such that a := >0- 


Thus we obtain 
|2 




a£iLMQ ^ A£2 


Il2(0,l) 


(28) 


a a 

and therefore the L^-norm of z is bounded by a constant. 

Now, let us multiply the first line of ( |22) l by xz and then 
integrate between 0 and L to get 

pL pL pL pL 

/ xzz'" dx + / xzz'dx + X xz^dx = / xzgdx 
Jo Jo Jo Jo 

After some integrations by parts, we get 


xzz'''dx = — 

3, 


zz''dx — 


xz'z''dx 


and 


— 2 ll^1lt=’(o,t) 


z'dx = —-I 


\l^(o,l)- 


Thus, plugging (30i and ([3T]l in (291, we obtain 


A\^'\\\'^(0,L) — c)\\A\'l'^{0,L) ^ 


xz^dx■ 


(30) 


(31) 


2 

L^{0,L) 


-2 

<M 


T 


xzg 


MI^(0,L) 


(32) 


where M is a constant which depends only on uo, L and a. 
In this way we see that 11-^^1112(0 is also bounded. 

From the Poincare inequality, we have the equivalence 
between ||-z'|li 2 (o 2 ,) ||z||^i(q^)- Thus, for any A > 0, 

a > 0 and y S L^{0,L), there exists M > 0 such that 
we have z G C, i.e. T{L^{Q, L)) C C. Then we apply the 
Schauder theorem (see e.g. [6, Theorem B.19]). Hence it 
concludes the proof of Lemma □ 

Thus, from [14, Theorem 4.2], using the Lemma [T] and 
A generates a semigroup of contraction T(t). With [1, 
Theorem 3.1], the proof of Theorem [T] is achieved. 

B. Asymptotic stability (Theorem 

The Lyapunov function related to (|^, which we will 
denote by E, is given by 

■E := ^ / y{t,xf‘dx, (33) 

and its derivative with respect to the time variable gives 

E < —a f y{t,x)sat{y{t,x))dx. (34) 

Jo 


y{t,x)sat{y{t,x)) >0 (35) 

Therefore we get 

^ < 0, (36) 

which means that, for all initial conditions in D{A), the 
solutions of (0 are stable. The attractivity has to be inspected 
too in order to finish the proof of the stability. 

Since we are in an infinite dimensional context, using the 
LaSalle’s Invariance Principle needs us to check whether the 
trajectories are compact. This precompactness is a corollary 
of the following lemma (which is very similar to [9, Lemma 
2], where a wave equation is considered). 

Lemma 4. The canonical embedding from D(A), equipped 
with the graph norm, into Lf(Q,L) is compact. 


Proof. Before proving this lemma, recall that its statement 
is equivalent to prove, for each sequence in D{A), which is 
bounded with the graph norm, that it exists a subsequence 
that (strongly) converges in L^{0,L). 

Let us recall the definition of the graph norm 



(37) 


Since for all (s, s) G C^, |s + sp < 2|sp + 2|sp, we get 
the following two inequalities 

||w||^(^) > 11^11^2(0^2.) (33) 

and 

r 

a , 


\d(A) > min I 1’ “ 


^ y I — asat{u)\‘^dx 
+ min y, J \u'" + u'+ asat{u)\‘^dx 


- ^ ^ ^ ^ |m'" + u'pdx. 


(39) 


Noticing that |lu"'||2= \\u"' + u' - we 

have 


h"'lli=(0.L) < 2|h"' + u'||i2(0,^) +2h'||i.(0_2,), (40) 


and using 

that = 

\W + W" 

- u'" + XU - 

^^Whio,Ly 

we obtain 




\\u\\Uo,L) 

<2||u 

+ L2(0,L) 

+ 

2||u"'- 

XU + XU\\\2(Q,L) 


<2||u' 

+ '^"'\\l2{0,L) 

+ 

4||u"'- 

^'‘^\\‘l^{0,L) 


+ 4| 

L2(0.L) 





<2||u' 

+ '^'"\\l2(0,L) 

+ 

^\W"\\Iho,l) 


L 


xu"'udx • 


■8||a;M||i2(0_i). 





From ( [30| l, we get 

rL 


'udx = -1 


Il 2 ( 0 ,l) 


and therefore 

lk'llL(0,L) <2h' + + 4h'"||i.(0.L) 


~ 12||u'|||2(o,l) 


■ SIIxmI 


L2(o,L)- 


Thus: 


(41) 


13|h'lli.(0.L) <2|k' + «'"lli3(0.L) + 4|k'"||i.(0.L) 

+ 8T^||u|||2(o,l) 

Plugging inequality (|40)i in (|42]i, we have 


(42) 


13|h'lli.(0,L) <2h' + «"'||i.(0,i) 

+ 4 {^\\u"' + M^||i2(o,L) + 2||M'|||2(o^i)^ 
+ 3F'^||u|li2(o_i) 

<10||tt' + «'"|li2(o.r.)+8||u'||i2(o.i) 

+ 8T^||u|||2(o,l)- 

and therefore 

0 7-2 

lh'llL(o,L) <2|k' + «'"||i2(o.i) + —||u||i2(o.,.) (43) 

Considering Equations ( [38l l and ( [39] l, it leads us to the 
following inequality, for all u G D{A) 


u'\\Uo,L) < A||urz,(^) 


(44) 


where A is a term which depends on L and a. 

Thus, if we consider now a sequence {u„}„gN in ^*(24) 
bounded for the graph norm of D{A), we have from (44 1 that 
this sequence is bounded in Hq{0,L). Since the canonical 
embedding from Hq{0,L) to L^(0, L) is compact, there 
exists a subsequence still denoted {«„}„£« such that Un ^ u 
in L^(0,L). Thus u belongs to L^{0,L), which concludes 
the lemma. □ 


Now we apply the LaSalle’s Invariance Principle. 

Using the fact that A generates a semi-group of contrac¬ 
tion, then from [1, Theoreme 3.1, Page 54], we get, for all 
t > 0 and for all ?/(0, •) € D{A), 


\\yi.t,-)\\L^O,L) < l|y(0,-)llL2(O.L) (45) 

and 

\\Ay{t, ■)\\l^(^o,l) < \\Ay{0,-)\\L^(o,L)- (46) 

Therefore, thanks to Lemma we see that the trajectory 
{v{t) = S{t)vo, f > 0} is precompact in L^{0,L), then the 
w-limit set •))] C D{A), is not empty and invariant 

to the nonlinear semigroup S{t) (see [22, Theorem 3.1]). 

Let us consider a strong solution such that E{t) = 0, for 
all f > 0. It follows from ( [3^ that y{t,x) = 0 for almost 
X in (0,L). Therefore the convergence property ( [^ h olds 
along the strong solutions to the nonlinear equation^. 


Using the density of D{A) and the existence of weak 
solutions, we end the proof by extending the result to any 
initial condition in L‘^{0,L). 

IV. Simulation 

Let us discretize the PDE (0 by means of finite difference 
method (see e.g. [21] for an introduction on the numerical 
scheme of a generalized Korteweg-de Vries equation). The 
time and the space steps are chosen such that the stability 
condition of the numerical scheme is satished. 

We choose L = 27r, yo{x) = 100(1 — cos(a;)) for all 
X G [0,27r] and a = 1. Let us numerically compute the 
solution of 0- On Figure there is no saturation in the 
dynamics. On Figure 1^ there is a saturation with a level 
Mo = 1. On Figure l^the feedback law is saturated with 
a level uq = 3. Figures and illustrate the evolution of 
the Lyapunov function E with respect to the time (without 
saturation and with a saturation level equals to 3). 
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Fig. 1. Solution y(t,x) with a feedback law without saturation 
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Fig. 2. Solution y(t,x) with a saturated feedback law and uq = 1. 
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Fig. 3. Solution y(t,x) with a saturated feedback law and no = 3. 




Fig. 4. Time evolution of the Fig. 5. Time evolution of the Lya- 

Lyapunov function E without sat- punov function E with a saturation 

uration. no = 3. 


V. Conclusion 

In this paper, we have studied the well-posedness and the 
asymptotic stability of a linear Korteweg-de Vries equation 
with a saturated distributed control. The well-posedness issue 
has been tackled by using the nonlinear semigroup theory 
and we proved the stability by using a sector condition 
and Lyapunov theory for infinite dimensional system. We 
illustrate our result on some simulations, which show that the 
smaller is the saturation level, the slower is the convergence 
to zero. 

To conclude, let us state some questions arising in this 
context: 

1. Can we extend our theorems to the nonlinear Korteweg-de 
Vries equation? 

2. As mentioned in the introduction, even if the internal 
control (without any constraints) acts only on a part of the 
domain, the stability still holds. Is it true with a saturated 
control? 

3. Can we recover the exponential stability with the saturated 
input? 

4. Can we apply the same method for other partial differ¬ 
ential equations? An interesting model could be the one¬ 
dimensional Kuramoto-Sivashinky equation. 
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